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ABSTRACT 

An analytically simple and tractable formula for the start-up autocovariances of periodic 
ARM A (PARMA) models is provided. 
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1. INTRODUCTION 

Autocovariance calculation procedures for PARMA models are generally carried out us- 
ing the periodic Yule Walker equations (e.g. Bentarzi and Aknouche, 2005). This approach 
has been considered earlier by Li and Hui (1988) for calculating PARMA autocovariances, 
where the (p+l) -start-up autocovariances, for all seasons, were given through a matrix equa- 
tion A'j = y, which is solved for 7 (7 being the (p + 1)S x 1-vector of the start-up autocovari- 
ances). The latter equation is however analytically and computationally intractable since 
the matrix A is not given explicitly but formed through an appropriate algorithm. Adopting 
the same approach, Shao and Lund (2004) showed that the r-start-up (r = max(p, q) + 1) 
autocovariances may be obtained by solving a linear system TUj = « for 7, where T and 
U are matrices of dimensions rS x (p + l)rS and (p + l)rS x rS, respectively. While these 
matrices are given explicitly, the method remains relatively cumbersome since it requires an 
increasing bookkeeping due to the matrix product. This note proposes an improved compu- 
tation procedure for calculating the PARMA autocovariances based on the latter approach. 
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The proposed method computes the (p + l)-start-up autocovariances based on a linear sys- 
tem with a corresponding matrix given explicitly, whose analytical form exhibits a circular 
structure, naturally assorted with the model periodicity. 



2. THE METHOD 

Consider a causal PARMA model of orders (p, q) and period S 

v <? 
yiffiyv-j+ns = y^B^Ey-j+nS, 1 < V < S, n G Z, (1) 

3=0 3=0 

where ^ = 6 { v) = -1 and {e t ,t G Z} is a periodic white noise process, i.e., a sequence of 
uncorrelated random variables with mean zero and variance E(e% +nS ) = a%, for 1 < v < S 
and jiGZ. 

Let 7^ = E {Uv+nsUv+nS-h) be the autocovariance function at season v and lag h G Z. 
Then, it is well known (Li and Hui, 1988; Shao and Lund, 2004) that multiplying (1) by 
y v+n s-h an d tacking expectation, the (lh^ are completely identified from the difference 
equation 

7?° - jyf^-3 j) = -E^Vjr fc fc) ^.i [fc < ff] , h > o, ( 2 ) 

3=1 j=h 

where the normalized cross- autocovariances (tp^), coefficients of the unique causal represen- 
tation of the PARMA process {yt,t G Z}, are given by (see e.g. Lund and Basawa, 2000; 
Shao and Lund, 2004) 

min(fc,p) 

V£° = -^%<«] + E ^t/\ k >1, v = 1, 5, (3) 

3=1 

with t^' 1 = 1 stands for the indicator function). 

Equation (2) needs to be started from the knowledge of 7^, < h < p and 1 < v < S. 
Once these start-up values are given, the 7^ for h > p may be obtained recursively from 
(2) while invoking (3). For 7^ with negative lags, we may use the well-known relation 
7^2 — ■ The mam result of this note is to formulate a linear system for computing the 

p + 1 necessary starting autocovariances. Let 7 = (7^, 7p*\ 7o 2 \ • 7p 2 \ 1q S \ Tp^)' 
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be the S(p + l)-vector of such values and £ be the S(p + l)-vector whose entries ChS+v = 
i2Qj V) i>\ V -h )a v-j, for 1 < v < s and < h < p, are the right-hand sides of (2). 

j=h 

Define the (p+ l)-square matrices ip^ (h > 0, 1 < v < S) and (k, v G {l,...,^}) 



as follows 
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and the S(p + l)-square matrix 



(i) 



(i) 



(2) 



5-1 
(2) 



(2) 



(2) 



for h = 0, ...,p 
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where mxn denotes the null matrix of dimension m x n. Then, the starting autocovariance 
vector 7 is the unique solution of the following linear system 



$7 = C, 



(6) 



whenever model (1) is causal. Note that, in view of (4), the infinite sum in (5) contains only 
p non zero terms. It may be possible to reduce the complexity of forming <fr using its circular 



property. Indeed, equation (5) may be used to only evaluate the first bloc 3>j^ , k 



.,S. 



The blocs <&j^ (v — 2, S) would be deduced from by substituting the corresponding 
parameters (fi^p by (fif* for j = I,..., p. 



(i) 



3. CONCLUDING REMARKS 

Despite the simplicity of the proposed method, it has the drawbacks that the starting au- 
tocovariances for all seasons are computed in the same bloc, thereby requiring O ((S(p + l)) 3 ) 
operations, which might be very costly for models with a large period. This is the main lim- 
itation of the periodic Yule Walker approach compared to which the methods that compute 
autocovariances for distinct seasons separately (e.g. Aknouche, 2007) are more suitable. 
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